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Abstract 

CN| ■ We present a mathematically rigorous analysis of the ground state of a dilute, interacting 

Bose gas in a three-dimensional trap that is strongly confining in one direction so that 
the system becomes effectively two-dimensional. The parameters involved are the particle 
number, N 1, the two-dimensional extension, L, of the gas cloud in the trap, the thickness, 
h <C L of the trap, and the scattering length a of the interaction potential. Our analysis 
starts from the full many-body Hamiltonian with an interaction potential that is assumed to 
be repulsive, radially symmetric and of short range, but otherwise arbitrary. In particular, 
hard cores are allowed. Under the premises that the confining energy, ~ l//i 2 , is much 
. larger than the internal energy per particle, and a/h — > 0, we prove that the system can be 

Oh' treated as a gas of two-dimensional bosons with scattering length a 2 D = h exp(— (const. )h/a). 

<-j ■ In the parameter region where a/h <C |ln(/3/i 2 )| , with p ~ N/L 2 the mean density, the 

system is described by a two-dimensional Gross-Pitaevskii density functional with coupling 
parameter ~ Na/h. If |ln(p/i 2 )| _1 < a/h the coupling parameter is ~ TVj hi(p/i 2 )| _1 and 
thus independent of a. In both cases Bose-Einstein condensation in the ground state holds, 
provided the coupling parameter stays bounded. 



X. 

1 Introduction 



In recent experiments dilute Bose gases have been confined in magneto-optical traps in such a 
way that the particle motion is essentially frozen in one or two directions and the system becomes 
effectively lower dimensional [8, 9, 33, 20, 23, 10, 31]. This is an intrinsically quantum mechanical 
phenomenon because it is not necessary to have a trap width or thickness that is the size of an 
atom, but it suffices that the energy gap for the motion in the strongly confined direction(s) 
is large compared to the internal energy per particle. The case of highly elongated ("cigar- 
shaped") traps has received particular attention [21, 4, 2, 3, 7, 26, 15] (see also [15] for further 
references), because it opens the possibility to realize the one-dimensional Lieb-Liniger model 
[11], and even the limiting Girardeau- Tonks case [6, 10, 23] that exhibits strong correlations. A 
detailed, rigorous derivation of the one-dimensional behavior from the many-body Hamiltonian 
of a three-dimensional gas was given in the paper [15]. This is not a simple problem, one reason 
being that an approximate factorization of the ground state wave function in the longitudinal 
and transverse variables is in general not possible (in particular not for hard core potentials), 
and the proofs in [15] are, in fact, quite long. 
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In the present paper we carry out a corresponding analysis for thin, disc-shaped traps, i.e., 
traps with strong confinement in one direction so that a two-dimensional behavior is expected. 
Experimental realizations of such systems and possible mechanisms for creating them are dis- 
cussed, e.g., in [8, 9, 1, 38, 31]. On the theoretical side the references [25, 22, 27, 28, 24] contain 
many valuable insights into their properties. There are several similarities with the emergence 
of one-dimensional behavior in cigar-shaped traps, but also some notable differences. Like for 
cigar-shaped traps, there is a basic division of the parameter domain into two regions: one where 
a limit of a three-dimensional Gross-Pitaevskii (GP) theory applies, and a complementary re- 
gion described by a "truly" low dimensional theory. In the case discussed in [15] the latter is a 
density functional theory based on the exact Lieb-Liniger solution for the energy of a strongly 
interacting (and highly correlated) one-dimensional gas with delta interactions. (Note that in 
ID strong interactions means low density.) In the present case, on the other hand, the gas is 
weakly interacting in all parameter regions. In the region not accessible from 3D GP theory the 
energy formula [32, 19, 14] for a dilute two-dimensional Bose gas with a logarithmic dependence 
on the density applies. To enter this region extreme dilution is required. The Lieb-Liniger region 
in the ID case demands also quite dilute systems, but the requirement is even more stringent in 
2D. This will be explained further below. 

We recall from [32, 19] that the energy per particle of a dilute, homogeneous, two-dimensional 
Bose gas with density /?2D and scattering length a2D of the interaction potential is (in units such 
that h = 2m = 1) 

e 2 D ~ 47rp2D| ln(p 2 Dal D )| _1 . (1.1) 
The corresponding result in three dimensions is 

e 3D « 47r/93 D a 3 D, (1.2) 

cf. [18]. In the following we shall denote the two-dimensional density, p2D, simply by p and the 
three dimensional scattering length, 030, by a. 

Our results establish rigorously, and in the many-body context, a relation between a, the 
thickness, h, of the trap (assumed to be 3> a) and the effective two-dimensional scattering length, 
a2D- Essentially, as h tends to zero, a2D = hexp(— (const. )h/a). (The precise formula is given in 
Eq. (1.17) below.) If \ln ph 2 \ <C h/a, then |ln(po| D )| « h/a, and the two-dimensional formula 
(1.1) leads to the same result as the three dimensional formula (1.2), because p^n ~ p/h. The 
"true" two dimensional region requires |ln/)/i 2 | > h/a and hence the condition p -1 / 2 > he h ^ a 
for the interparticle distance, p~ 1 / 2 . This should be compared with the corresponding condition 
for the ID Lieb-Liniger region in [15] where the interparticle distance is "only" required to be 
of the order or larger than h 2 /a. 

The basic formula a2D = hexp(— (const. )h/a) for the scattering length appeared, to the best 
of our knowledge, first in [25]. It can be motivated by considering a weak, bounded potential, 
where perturbation theory can be used to compute the energy for a two-body problem that is 
directly related to the scattering length, cf. Appendix A in [19]. This perturbative calculation 
is carried out in Section 4 as a step in the proof of a lower bound for the many-body energy; 
its relation to the formula for a2D is explained in the Remark after Corollary 4.2. We wish to 
stress, however, that deriving this formula in the context of two-body scattering is only a step 
towards the solution of the many-body problem that is the concern of the present paper. 

We now define the setting and state the results more precisely. We consider N identical, 
spinless bosons in a confining, three-dimensional trap potential and with a repulsive, rotationally 
symmetric pair interaction. We take the direction of strong confinement as the z-direction and 
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write the points x £ I 3 as (x, z), x £ 1R 2 , z£R. The Hamiltonian is 

N 

H N , L , h , a = J2(- A i + V L,h(Xi))+ "aCIXi-Xj-D (1.3) 

i=l l<i<j<N 

with 

Vb lfc (x) = ^) + y/w = ^( rt )4 yl(rlz) ' (L4) 

«a(|x|) = ^(a^lxl). (1.5) 

The confining potentials V and y 1 - are assumed to be locally bounded and tend to oo as \x\ 
and \z\ tend to oo. The interaction potential v is assumed to be nonnegative, of finite range 
and with scattering length 1; the scaled potential v a then has scattering length a. We regard v, 
V 1 - and V as fixed and L,h,a as scaling parameters. The Hamiltonian (1.3) acts on symmetric 
wave functions in L 2 (M. 3N , dxi ■ ■ ■ dxjy). Its ground state energy, E® M (N, L, h, a), scales with L 
as 

E QM (N, L, h, a) = -^E QM (N, 1, h/L, a/L). (1.6) 



Taking N — > oo but keeping h/L and Na/L fixed leads to a three dimensional Gross-Pitaevskii 
description of the ground state as proved in [13]. The corresponding energy functional is 

4dM= / {|V0(x)| 2 + ^(x)| < A(x)| 2 +4vriVa|0(x)| 4 }d 3 x (1.7) 

and the energy per particle is 

E^(N,L,h,a)/N = inf{£ 3 G D P M: / |0(x)| 2 d 2 x = 1} 

= (1/ L 2 )E^ (1,1, h/L, Na/L). (1.8) 

By Theorem 1.1 in [13], we have, for fixed h/L and Na/L, 

E^(N.L M = 
N^oo E^(N,L,h,a) V ; 

It is important to note, however, that the estimates in [13] are not uniform in the ratio h/L and 
the question what happens if h/L — > is not addressed in that paper. It will be shown in the 
next section that a part of the parameter range for thin traps can be treated by considering, at 
fixed Na/h, the h/L —> limit of £^(1, 1, h/L, Na/L), with the ground state energy for the 
transverse motion, ~ 1/h 2 , subtracted. But this limit can evidently never lead to a logarithmic 
dependence on the density and it does not give the correct limit formula for the energy in the 
whole parameter range. 

To cover all cases we have to consider a two-dimensional Gross-Pitaevskii theory of the type 
studied in [14], i.e., 

£ 2 dV]= / {\V^x)\ 2 + V L (x)\^x)\ 2 +A^Ng\ V {x)t}d 2 x (1.10) 

with 

g=\ln(pa 2 D \-\ (1.11) 
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Here p is the mean density, denned as in Eq. (1.6) in [14]. An explicit formula that is valid in 
the case Ng » 1 can be states as follows. Let 

p7(x) = ^[»7-Vl(x)} + (1.12) 

be the 'Thomas Fermi' density for N particles at coupling constant 1 in the potential Vl, where 
pj^ is chosen so that f pjf = N. Then 

p = N- 1 [ pl F (x) 2 dx. (1.13) 
Jm. 2 

For simplicity we shall assume that V is homogeneous of some degree p > 0, i.e., V(Xx) = 
X p V(x), and in this case 

p~ atp/(p+2)/ L 2 = N /l2 with i = N i/(p+2) L _ ( L14 ^ 

The length L measures the effective extension of the gas cloud of the N particles in the two- 
dimensional trap. A box potential corresponds to L = L, i.e., p = oo and hence p ~ N/L 2 . 

The case Ng = 0(1) requires a closer look at the definition of p. First, for any value of Ng 
we can consider the minimizer of (1.10) with normalization H2 = 1- The corresponding 
mean density is 

„GP|4 



PNg = N J \^ g \\ (1.15) 

A general definition of p amounts to solving the equation p = pNg with g as in (1.11). As 
discussed in [14] this gives the same result as (1.13) to leading order in g when Ng » 1. In the 
case I In Nh/L\ <C h/a (referred to as 'Region I' below) the coupling constant is simply a/h and 
thus independent of N. Moreover, in a homogeneous potential of degree p the effective length 
scale L is ~ (Ng + l) l /^ + ^L and thus of order L if Ng = 0(1). 
The energy per particle corresponding to (1.10) is 

E^(N,L,g)/N = inf{£ 2 G D P M : / \^{x)\ 2 d 2 x = 1} = (1/L 2 )£ 2 G D P (1, l,Ng). (1.16) 

Let Sh be the normalized ground state wave function of the one-particle Hamiltonian 
—d 2 /dz 2 + Vfj-(z). It can be written as Sh{z) = h~ 1 / 2 s(h~ 1 z) and the ground state energy 
as = /i _2 e _L , where s(z) and e 1 - are, respectively, the ground state wave function and ground 
state energy of —d 2 /dz 2 +V- L (z). We define the two dimensional scattering length by the formula 

a 2D = /iexp (-(/ s{z) 4 dzy 1 h/2a) . (1.17) 

Then, using (1.11), 

g=\- \n{ph 2 ) + (/ s(zfdz)- 1 h/a\- 1 . (1.18) 
The justification of the definition (1.17) is Theorem 1.1 below. 

Remark. Since a2D appears only under a logarithm, and a/h — > 0, one could, at least as far 
as leading order computations are concerned, equally well define the two dimensional scattering 
length as a 2D = b exp (—(J s(z) A dz)~ 1 h/2a\ with b satisfying ca < b < C h for some constants 

1 2d) 



c > 0, C < 00. In fact, if g' = \ ln(p(a' 2D ) 2 )\ 1 , then 



+ 21n ^) ,1 ( ll 9 ) 

g' I - ln(ph 2 ) + (const. )h/a\ 



because (a/h)ln(b/h) — > 0. 

We can now state the limit theorem for the ground state energy: 



5 



Theorem 1.1 (Prom 3D to 2D, ground state energy). Let N — > oo and at the same time 
h/L — > and a//i — > in snc/i a way i/iai /i 2 po — > (with g given by Eq. (1.18) ). Then 

E° P (AT,L, 9 ) 

Remarks: 1. The condition /i 2 /?<7 — > means that the ground state energy h~ 2 e^ associated 
with the confining potential in the z-direction is much larger than the energy pg. This is the 
condition of strong confinement in the z-direction. In the case that h/a^> \\\i(ph 2 )\ we have 
g ~ a/h and hence the condition in that region is equivalent to 

M«l. (1.21) 

On the other hand, if h/a < |ln(/)/i 2 )| the strong confinement condition is equivalent to 
h 2 p\ \n.{h 2 p)\~ l <C 1, which means simply that 

p7i 2 «l. (1.22) 

Both (1.21) and (1.22) clearly imply pa 2 D <C 1, i.e., the gas is dilute in the 2D sense (and 
also in the 3D sense, pzoo? <C 1, because p^D = p/h). This is different from the situation in 
cigar-shaped traps considered in [15] where the gas can be either dilute or dense in the ID sense, 
depending on the parameters (although it is always dilute in the 3D sense). 

2. It is, in fact, not necessary to demand h/L — > explicitly in Theorem 1.1. The reason is 
as follows. In the region where h/a < |ln(/?/i 2 )|, the strong confinement condition ph 2 <C 1 
immediately implies h/L <C 1 because p S> l/L 2 , cf. Eq. (1.14). If h/a 3> |ln(p7t 2 )|, then 
at least pah <C 1 holds true. This leaves only the alternatives h/L — > 0, or, if h/L stays 
bounded away from zero, Na/L — > 0. But the latter alternative means, by the three dimensional 
Gross-Pitaevskii limit theorem [13], that the energy converges to the energy of a noninteracting, 
trapped gas, for which (1.20) obviously holds true. 

We shall refer to the parameter region where h/a ^> |ln(p7i 2 )| as Region I, and the one 
where h/a < | ln(p^, 2 ) | as Region II. In Region I we can take 

g = (f s {zfdz)a/h. (1.23) 

In Region II g ~ | ln(p/i 2 )| _1 , and in the extreme case that h/a <C | ln(p/i 2 )|, 

g=\Hph 2 )\- 1 . (1.24) 

In particular g is then independent of a (but dependent on p). As remarked earlier, Region II 
applies only to very dilute gases since it requires interparticle distances p -1 / 2 > he h l a . 

By Eq. (1.16) the relevant coupling parameter is Ng rather than g itself, and both Region 
I and Region II can be divided further, according to Ng <^ 1, Ng ~ 1, or Ng 3> 1. The case 
Ng <C 1 corresponds simply to an ideal gas in the external trap potential. Note that this limit 
can both be reached from Region I by taking a/h — > at fixed ph 2 , or from Region II by letting 
ph 2 tend more rapidly to zero than e~ h l a . The case Ng ~ 1 in Region I corresponds to a GP 
theory with coupling parameter ~ Na/h as was already explained, in particular after Eq. (1.9). 
The case Ng 3> 1 is the 'Thomas-Fermi' case where the gradient term in the energy functional 
(1.10) can be ignored. In Region II, the case Ng < 1 requires p -1 / 2 > he N and is thus mainly 
of academic interest, while p -1 / 2 <C he N (but still he h l a < p^ 1 ^ 2 ) corresponds to the TF case. 
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The subdivision of the parameter range just described is somewhat different from the 
situation described in [15]. The reason is the different form of the energy per particle of the low 
dimensional gas as function of the density. 

Theorem 1.1 is a limit theorem for the ground state energy. By standard arguments (variation 
with respect to the trapping potential) it implies a convergence result for the one particle density 
in the ground state, PmlhJ^) Section 2 in [15] and Theorem 8.2 in [17]): 

Define the 2D QM density by integrating over the transverse variable z, i.e., 

pTl,h^)--= J PN M LAa( x >z)dz • (1-25) 

With L the extension of the system in the 2D trap, cf. (1.14), define the rescaled GP density p 
by 

p(x) = L 2 \<p GP (Lx)\ 2 (1.26) 

where ip GP is the minimizer of (1.10) with normalization f \tp GP (x)\ 2 dx = 1. (Note that p 
depends on N, L and g.) 

Theorem 1.2 (2D limit for the density). In the same limit as considered in Theorem 1.1, 

(^pT L , h , a {Lx) - /?(*)) =0 (1.27) 

in weak L 1 sense. 

In the GP case where the coupling parameter Ng stays bounded a much stronger result can 
be proved, namely convergence of the 1-particle density matrix and Bose-Einstein condensation 
(BEC) in the ground state. Recall that the one-body density matrix obtained from the ground 
state wave function is 

7tt (x,x') = N / ^o(x,X2,...,x A r)^o(x / ,X2,...,x A r)*dx2---dx A r. (1.28) 

JR3(AT-1) 

BEC means that in the ./V — ► oo limit it factorizes as Nip(x)ifj(x') for some normalized ip. This, 
in fact, is 100% condensation and was proved in [12] for a fixed trap potential in the Gross- 
Pitaevskii limit, i.e., for both h/L and Na/L fixed as N — > oo. Here we extend this result to the 
case h/L — > with Ng and L fixed. The function tp is the minimizer of the 2D GP functional 
(1.10) times the transverse function Sh(z): 

Theorem 1.3 (BEC in GP limit). // N — > oo, h/L — > while Ng and L are fixed, then 

-^7*o(*, hz; x', hz') -> <p GP {x)ip GP {x')s(z)s(z') (1.29) 

in trace norm. Here ip GP is the normalized minimizer of the GP functional (1.10). 

There is a variant of Theorem 1.1 that applies to the thermodynamic limit in the 2D variable 
x where the density becomes homogeneous in this variable. Let E^, ox (N,L,a) be the ground 
state energy of the Hamiltonian (1.3) with the potential Vl(x) replaced by a 2D box of side 
length L and define, for fixed p, h and a, 

, , , E hox (N,L,a)-Nh- 2 e ± 
e 2D (p,h,a) := hm . (1.30) 

N,L^oo,N/L 2 =p iV 
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This is the energy per particle in the 2D thermodynamic limit (with the confining energy sub- 
tracted) and by standard arguments (cf. e.g., [30]) it is independent of the boundary conditions 
(Dirichlet, Neumann or periodic) imposed in the definition of E^ ox (N, L, a). We then have 

Theorem 1.4 (From 3D to 2D, homogeneous case). If h 2 pg — > and a/h — > 0, with g 
given by Eq. (1.18), then 

lim e 2D (p,M) =1 (L31) 
iirpg 

The essentials of the dimensional reduction are already contained in the proof of Theorem 1.4 
and are more transparent in this case than for inhomogeneous gases in Theorem 1.1. Hence we 
shall focus on the proof of Theorem 1.4, while the modifications that have to be made for a proof 
of Theorem 1.1 will be more briefly described with appropriate references to [13, 14, 15, 36], 
where analogous problems are discussed in some detail. 

An abbreviated version of this paper appears as Chapter 9 in the Oberwolfach Seminars 
volume [17]. 



2 The 2D limit of 3D GP theory 

As in [15], certain aspects of the dimensional reduction of the many-body system can be seen 
already in the much simpler context of GP theory. We therefore begin by considering the 
h/L — > limit of the 3D GP ground state energy. The result is, apart from the confining energy, 
the 2D GP energy with coupling constant g ~ a/h. This shows in particular that Region II, 
where g ~ | ln(pft, 2 )| -1 , cannot be reached as a limit of 3D GP theory. 

Theorem 2.1 (2D limit of 3D GP energy). Define g = (/ s{z) A dz) a/h. If h/L -» 0, then 

E^(N,L,h,a)-Nh~ 2 e 1 - 

E%(N,L,g) {ZA) 

uniformly in the parameters, as long as pah — ► 0. 

Remark. Since E^ (1, L, Ng) ~ L~ 2 + pa/h, the condition pah — > is equivalent to 
h 2 E^(l,L,Ng) -» 0, which means simply that the 2D GP energy per particle is much less 
than the confining energy, ~ 1/h 2 . 

Proof. Because of the scaling relation (1.8) it suffices to consider the case N = 1 and L = 1. 
For an upper bound to the 3D GP ground state energy we make the ansatz 

<Hx) = <p GP (x)s h (z), (2.2) 

where ip GP is the minimizer of the 2D GP functional with coupling constant g. Then 

£ 3 G dM =e ± /h 2 + E<£(l,l,g) (2.3) 

and hence 

4 P (l,l ) ft,o)-e i // l 2 < £ 2 G D P (1,1, 5 ). (2.4) 
For the lower bound we consider the one-particle Hamiltonian (in 3D) 

H Ka = -A + ^(x) + SH^gpWI 2 ^) 2 . (2-5) 
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Taking the 3D GP minimizer $ as a test state gives 

infspec-Hha < 1, h, a) - 4wa / |$(x)| 4 d 3 x 

+8vra y |^Gp(x)| 2 s/ l (z) 2 |$(x)| 2 d 3 x 
< £ 3 G D P (1,1,M) +4vra / |y? GP (x)| 4 S/l (z) 4 d 3 x 
= Ejjg(l,l,h,a)+4wg ! \<p GP (x)\ 4 dx. (2.6) 

To bound a from below we consider first for fixed x£l 2 the Hamiltonian (in ID) 

H h , a , x = -d 2 z + v£(z) + 8na\ V GP {x)\ 2 s h {z) 2 . (2.7) 

We regard — d 2 + Vf^(z) as its "free" part and 8ira\(p GP (x)\ 2 Sh{z) 2 as a perturbation. Since the 
perturbation is positive all eigenvalues of H^x are & t least as large as those of — d 2 + Vh(z); in 
particular the first excited eigenvalue is ~ 1/ti 2 . The expectation value in the ground state Sh 
of the free part is 

(H h ,a,x) = e ± /h 2 + 87rg\cp GP (x)\ 2 . (2.8) 
Temple's inequality [37, 29] gives 

H lw > (e^ + ^ (,)P) (l - jg^A^> ) , ,2.0, 



where e^/h 2 is the lowest eigenvalue above the ground state energy of — <9^ + V^-(2:). Since 



(2.10) 



H Ka , x s h = (e ± /h 2 )s h + 87ra|99 GP (x)| 2 4 
we have (H hj<ljX - (H h>ajX ))s h = 8Tr\^ GP (x)\ 2 (as 3 h - gs h ) and hence, using g = a js\ = (a/h) f s 4 , 

{(Hh,a,x — (Hh,a,x)) 2 } = (8vr) 2 |^ GP (x) | 4 J (a3 h {z) 3 ~ 9 S h (z)) 2 dz 



< (8,) 2 y GP \t(a/h) 2 

< const. £ 2 GP (1,1,<?) 2 



(2.11) 



where we have used Lemma 2.1 in [14] to bound the term g\\ip H 2 ^ by const. £^(1, 1, g). We 
thus see from (2.8) and the assumption h 2 E G ^(l, l,g) ^ that the error term in the Temple 
inequality (2.9) is o(l). 

Now Hh a = — A x + V(x) + Hh a j., so from (2.9) we conclude that 



H h , a > ( (e 4 -/^) - A, + V(x) + 8irg\v GP (x)\ 2 ) (1 - o(l)). 



(2.12) 



On the other hand, the lowest energy of -A x + V(x) + 87rg\<p GP (x)\ 2 is just E GP (l,l,g) + 
Airg J R2 \<PGp(x)\ 4 dx. Combining (2.6) and (2.12) we thus get 



£ 3 G D P (1, 1, h, a) - e^/tf > E%(1, 1, g)(l - o(l)). 



.J- /i,2 



GP/ 



(2.13) 
□ 
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3 Upper bound 
3.1 Finite n bounds 

We now turn to the many-body problem, i.e., the proof of Theorems 1.1 and 1.4. Like in [15] 
the key lemmas are energy bounds in boxes with finite particle number. The bounds for the 
total system are obtained by distributing the particles optimally among the boxes. We start 
with the upper bound for the energy in a single box. 
Consider the Hamiltonian 

n 

# = E(- A * + ^))+ E "adXi-X,-!) (3.1) 

i=l l<i<j<n 

in a region A = A2 x R where A2 denotes a box of side length i in the 2D x variables. For the 
upper bound on the ground state energy of (3.1) we impose Dirichlet boundary conditions on 
the 2D Laplacian. The goal is to prove, for a given 2D density p and parameters a and h, that 
for a suitable choice of I and corresponding particle number n = p£ 2 the energy per particle, 
with the confining energy e 1 - /h 2 subtracted, is bounded above by 

47rp|ln(p4 D )|- 1 (l+o(l)) (3.2) 

where a2D is given by Eq. (1.17). Moreover, the Dirichlet localization energy per particle, ~ l/^ 2 , 
should be small compared to (3.2). The relative error, o(l), in (3.2) tends to zero with the small 
parameters a/h and pah (Region I), or a/h and ph 2 (Region II). 

The choice of variational functions depends on the parameter regions and we are first con- 
cerned with the Region II, i.e., the case | ln(ph 2 )\ > h/a. 

3.1.1 Upper bound in Region II 

Let fo(r) be the solution of the zero energy scattering equation 

-A/ + ±v a f = 0, (3.3) 

normalized so that fo(r) = (1— a/r) for r > Rq with Rq the range of v a . Note that Rq = (const.) a 
by the scaling (1.5). The function /o satisfies < fo(r) < 1 and < /o(r) < min{l/r, a/r 2 }. 
This is seen by writing fo(r) = u(r)/r and fo(r) = v! \r) / V — u(r) / 'r 2 with u"{r) = ^v(r)u(r) > 0. 
Since u(r) = r — a for r > Rq and u(0) = 0, convexity implies u(r) > max{0, r — a} and u'(r) < 1. 
Hence < /o(r) < min{l/r, a/r 2 } and < fo(r) < lim r ^oo(l — a/r) = 1. 
For R > Rq we define 

f(r) = / (r)/(l - a/R) for < r < R, and f(r) = 1 for r > R. (3.4) 

We define a two-dimensional potential by 

W {x) = 3M f [ /( | X |)2 + i^dxD/dxl) 2 ] dz. (3.5) 

Clearly, W(x) > 0, and W is rotationally symmetric with W(x) = for |x| > R. Moreover, by 
partial integration, using (3.3), it follows that W 6 L X (R 2 ) with 

/ W(x)dx= 8 ™l S ^ (l-a/R)-\ (3.6) 
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Define, for b > R, 

' ln(R/a 2D )/ln(fe/a 2D ) if < r < R 

ip{r) = I ln(r/a 2 D)/ln(6/a 2 D) if R<r <b 

1 if b<r 

As test function for the three dimensional Hamiltonian (3.1) we shall take 

*(xi,...,x n ) = F(xi,...,x n )G(xi,...,x n ) 

with 



(3-7) 



(3.8) 



F(xi, . . . ,x n ) = JJ/flxj - Xj|) and G(xi, . . . ,x n ) = JJ ip(\xi - xj\) JJ s h {z k ). (3.9) 

i<j i<j k=l 

The parameters R, b and also £ will eventually be chosen so that the errors compared to the 
expected leading term in the energy are small. 

As it stands, the function (3.8) does not satisfy Dirichlet boundary conditions but this can 
be taken care of by multiplying the function with additional factors at energy cost ~ 1/i 2 per 
particle, that will turn out to be small compared to the energy of (3.8). 

Since /(| |)(^(|xj - Xj\) = 1 for \x,i — Xj\ >b and s h is normalized, the norm of V can 

be estimated as in Eq. (3.9) in [15]. 



2n 



im{n — 1) b 2 
2 P 



(3.10) 



Next we consider the expectation value of H with the wave function ty. By partial integration 
we have, for every j, 

J |V,(FG)| 2 

= J G 2 \VjF\ 2 - J F 2 GAjG = J G 2 \VjF\ 2 - J F 2 Gd 2 Zj G- J F 2 G(aJg) 

= J G 2 \VjF\ 2 - J F 2 Gd 2 Zj G + J F 2 |vjGj 2 + 2 J FG(vJf) ■ (vJg) (3.11) 

where A!j and v!j are, respectively, the two dimensional Laplace operator and gradient. The 
term — f F 2 Gd 2 .G together with jV^F 2 G 2 gives the confinement energy, (e- L //t 2 )||vD r || 2 . 

Next we consider the first and the third term in (3.11). Since < / < 1, /' > and is 
normalized, we have 

£/f 2 |vJ!g| 2 + £ J \V jF \ 2 G 2 <Y: J |v5<I>| 2 + 2^ / '(ix,-x,|) 2 G 2 

3 3 3 i<3 

+ 4 E / / / (|x fc -x,|)/'(|x fe -x j |)G 2 



(3-12) 



k<i<j ' 



where we have denoted Yii<j <f(\ x i ~ x j\) by & for short. Moreover, since < ip < 1, 



2 f'(\*i ~ *i \?G 2 < 2 f'(\*i ~ X J \) 2 Sh(*i) 2 *i (*3 



(3.13) 



i<3 



i<j 
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By Young's inequality 

|4 



2/ fX\x i -y Lj \f Sh {z i ) 2 s j {z j ) 2 dz i dz j <' 2 ^! f'{\{Xi-Xj,z)\) 2 dz. (3.14) 

The right side gives rise to the first of the two terms in the formula (3.5) for the two dimensional 
potential W. The other term is provided by J F 2 G 2 v a {^-i — Xj), using < / < 1, < </? < 1 and 
Young's inequality. 

Altogether we obtain 

(*|fT|*> -{ne^/h 2 )^) <V I |v"$| 2 + V/ W(x i -x j )^ 2 + TZ 1 + TZ 2 (3.15) 

J I A n ' / A Tl 



with 



and 



K! = 2 f FG(vjF) • (VJG) (3.16) 
^2 = 4 V / / , (|x fc -x l |)/ / (|x fe -x,|)G 2 (3.17) 

k<i<j JAn 

< | n (n _ i)( n _ 2)^-3) / /'(| Xl - x 2 |)/'(|x 2 - M>h(zifs h {z 2 ) 2 s h {z,) 2 , 

where < if < 1 has been used for the last inequality. 

The error term TZ\ is easily dealt with: It is zero because ip(r) is constant for r < R and 
f(r) is constant for r > R. 

The other term, 7^-2 , is estimated as follows. Since f'(r) = for r > R we can use the Cauchy 
Schwarz inequality for the integration over xi at fixed x 2 to obtain 



/'(|xi -x 2 \)8 h (z 1 ) 2 dx 1 



/ 

< f / /'(| Xl -X 2 |) 2 d Xi y 2 ( / S h (zi) 4 d^) 
\J J \^|xi-x 2 |<H / 

< (47!-||s|| 00 a / J R 3 /3/i 2 ) 1/2 (3.18) 

with a' = a(l — a/R)^ 1 . The same estimate for the integration over X3 and a subsequent 
integration over x 2 gives 

1Z 2 < (constO^n 3 ^. (3.19) 

We need TZ 2 /(^>\^>) to be small compared to the leading term in the energy, ~ n 2 ^~ 2 | ln(/9/j 2 )| _1 
with p = n/l 2 . (Recall that we are in Region II where | ln(p/i 2 )| > h/a.) Moreover, the leading 
term should be large compared to the Dirichlet localization energy, which is ~ n/l 2 . We are 
thus lead to the conditions (the first comes from (3.10)): 

n 2 b 2 na'R 3 1 ln(p/i 2 ) I n 

— «i. — n^* 1 . ( 3 - 2 °) 

which can also be written 
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These conditions are fulfilled if we choose 

R = h, b = p- l / 2 \\n(ph 2 )\- a (3.22) 

with a > 1/2 and 

p- 1 / 2 ] ln^/i 2 )! 1 / 2 < £ < p- 1 / 2 ] Mph 2 )\ a . (3.23) 
Note also that n = pi 2 > 1. 



It remains to compare 



| V I!<D| 2 + W WO*-*,-)* 2 ] (3.24) 

with the expected leading term of the energy, i.e., A7r(n 2 /£ 2 )\ ln(na 2 D /.£ 2 )| _1 . 
We consider first the simplest case, i.e., n = 2. We have 

/ \V^\ 2 = (ln(b/a 2D )y 2 27r f b — = (\n(b/a 2 D)y 2 2TTln(b/R), (3.25) 

Inserting the formula (1.17) for a2D and using R = h, b = p~ l / 2 \\n{ph 2 )\- a and a' = a(l + o(l)) 
we have 

/ ( |vVi 2 + i^ 2 ) = 

27r(ln(6/a 2D ))~ 2 [ln(6//i) + (fe/2o'||s||t)] = 4vr| ln^alo)^^ + o(l)). (3.27) 
For n > 2 we can use the symmetry of to write, using (3.27) as well as < ip(r) < 1, 

£ / i v i l(|, i 2 + E / w(^-^ 2 

= nl f |vf<I>| 2 + ±V/ W^-X!)* 2 ] 

< 4Trn 2 £ 2( - n -^ \ \n(nal D /£ 2 )\-\l + o(l)) + ^ 3 (3.28) 

with 

7l 3 = n 3 ^- 3 ) / ^'(Ixa-xilVdara-xil). (3.29) 

We estimate TZs in the same way as (3.18), obtaining 

K 3 < (const.)£ 2( - n - 2 ^n 3 b 2 (\n(b/a 2 D)y 2 2TTln{b/R). (3.30) 

The condition that IZ3 has to be much smaller than the leading term, given by 
47m 2 ^ 2 ( n - 1 )|ln(na 2 , D /^ 2 )|- 1 , is equivalent to 

^-In(6/ J R)<1. (3.31) 
With the choice (3.22) this holds if a > 1/2. 
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3.1.2 Upper bound in Region I 

In Region I the ansatz (3.8) can still be used, but this time we take b = R, i.e., ip = 1. In 
this region (a/h) \ ln(ph 2 )\ = o(l) and the leading term in the energy is ~ n 2 £~ 2 a/h. Conditions 
(3.20) are now replaced by 

n 2 R 2 _ nR? na 

— <K1 - W* 1 ' ~h * (332) 

where have here used that a' = a(l + o(l)), provided R 2> a. Note that the last condition in 
(3.32) means in particular that u>l. Putting again p = n/£ 2 , (3.32) can be written as 

pH 2 R 2 ^l, ^-«1, (3.33) 
h pi z a 

By assumption, a/h <C 1, but also pah <C 1 by the condition of strong confinement, c.f. (1.21). 
We take 

R = a(pah)' 13 (3.34) 

with < (3, so R » a. Further restrictions come from the conditions (3.33): The first and the 
last of these conditions imply together that 

- < p£ 2 < (3.35) 
a pR 2 

which can be fulfilled if 

pah < {pah) 213 . (3.36) 
i.e., if (3 < \. Note that this implies in particular R p~ x l 2 . We can then take 

^ = p- 1 /2( /l / a )l/2 (pa/l) - 7 (3.37) 

with 

< 7 < (3.38) 
The second of the conditions (3.33) requires that 

^- = (pah) (a/h) 2 (pah)- 313 < 1, (3.39) 
which holds in any case if f3 < 1/3. A possible choice satisfying all conditions is 

"=5- ~< = T2- < 3 ' 40 > 

The error terms (3.33) are then bounded by (pah) 1 / 6 (first and third term) and (a/h) 2 (second 
term). 

Finally, with $ = 1, Eqs. (3.15), (3.10) and (3.6) give 

(^r^liTI*) - (ne^/h 2 ) < a|l ^ (l + o(l)). (3.41) 
This completes the proof of the upper bound in boxes with finite n. 
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3.2 Global bound, uniform case 

The upper bound for the energy per particle in the 2D thermodynamic limit needed for Theorem 
1.4 is now obtained by dividing R 2 into Dirichlet boxes with side length t satisfying (3.23) in 
Region II, or (3.35) in Region I, and distributing the particles evenly among the boxes. In other 
words, the trial wave function in a large box of side length L is 

a 

where a labels the boxes of side length i contained in the large box, and ^> a is the Dirichlet 
ground state wave function for n = pi 2 particles in the box a, with p = N/L 2 . The choice of I 
guarantees in particular that the error associated with the Dirichlet localization is negligible. In 
order to avoid contributions from the interaction between particles in different boxes the boxes 
should be separated by the range, i?o of the interaction potential and in the "corridors" between 
the boxes the wave function is put equal to zero. The number of particles in each box is then 
not exactly p£ 2 , but rather the smallest integer larger than pt 2 (lj(l + i?o)) 2 . In order that the 
errors are negligible one needs Rq/1 = o(l), as well as p£ 2 S> 1, and both are garanteed by the 
choice (3.37) of I. 

3.3 Global bound in a trap 

In Theorem 1.1 the system is inhomogeneous in 2D because of the trapping potential Vl(x), 
and the distribution of the particles among the boxes has to be adjusted to the density given by 
the GP minimizer (p GP (x). As long as £, given by (3.37) or (3.23), is small compared to L this 
can be done in an analogous way as in [15], Section 4.2. Since, however, the condition I < L 
requires Na/h » (pah)~ l l 12 in Region I, and N/\ ln(ph 2 )\ 3> |ln(p/i 2 )| e in Region II, this will 
not work in the whole parameter range. 

A better and only slightly more complicated choice that works in all cases is to replace the 
Jastrow-type function IL<j /(l x « ~~ x jl) by a Dyson wave function of the type considered in [5] 
and [13], i.e., 

JV 

F(xi,...,xjv) = J]/(t i (xi,...,x i )) (3.43) 

i=i 

where 

i»(xi, . . . ,Xj) = min |xj-x,-| (3.44) 

l<j<i— 1 

is the distance between Xj and its nearest neighbor among the points xi,...,Xj_i. To take 
care of the inhomogeneity in the 2D variables a factor rii l / 9GP ( x j) i s included, where (p GP the 
normalized minimizer of the 2D GP functional with coupling constant iVg. Moreover, the local 
behavior in the 2D variables is modelled by a Dyson wave function 

$(xi, . . . , x N ) = if(t' i (x 1 , x N )) (3.45) 

with 

t'i(xi, . . . , Xi) = min \xi — Xj\. (3.46) 
l<j<i— 1 

The global trial wave function is thus of the form 



*(xi, . . . ,Xjv) = F(xi, . . . ,Xjv)G(xi,. . . ,xjv) 



(3.47) 
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with F given by (3.43), and 

G(xi, . . . ,xjv) = $(xi, . . . , xjv) Yl V GP (xi)s h (zi). (3.48) 

The advantage of the Dyson wave functions over the Jastrow-type functions is that when comput- 
ing the expectation values (^f,H^f)/(^f,^) cancellations between numerator and denominator 
take place that effectively improve the estimate (3.10) of the norm of *S>. 



3.3.1 Region I 

In Region I, where one can take ip = 1, the computation of the expectation value of the Hamil- 
tonian is the same as in [13]. In fact, there is even no need to do the computations explicitly 
because the required bound can be obtained by combining Theorem 2.1 with Theorem III. 1 in 
[13]. One checks that the parameter b in Eq. (3.42) in [13] is proportional to p -1 / 3 /?, 1 / 3 and 
hence the error terms in the upper bound of E® M in terms of the 3D GP energy are of the 
order a/b = (a//i) 2 / 3 (pa/i) 1 / 3 . On the other hand, by Theorem 2.1 the 3D GP energy minus 
confining energy approaches, as h/L — > 0, the 2D GP energy with coupling constant (J s )a/h 
and the approximation is uniform in the parameters as long as {pah) — > 0. Altogether we have 
the bound 

E QM {N,L,h,a) - Nh~ 2 e ± < E^(N,L,{J s 4 )a/h){l + o{l)) (3.49) 
where the error term o(l) tend to zero if h/L — > 0, a/h — > and pah — > 0. 



3.3.2 Region II 

In Region II one uses the ansatz (3.47)— (3.48) with / as in (3.4) with R = h and ip as in (3.7) 
with b = /? _1 / 2 | ln(p/i 2 )[ a . Since the computations are analogous to those in [13] and in Eqs. 
(3.11)-(3.29) above we shall not write all details explicitly. Proceeding as in Section 3.1 in [13], 
using (3.11), we obtain for the expectation value of the Hamiltonian (1.3) 

(*|*) - 2^ J\^\ 2 ^ Jl^l 2 

+2 E (350) 

k<i<j ■' ' 

" SF 2 G(-V 2 G) + \^V L , h (^) 
1 for i = k 

— 1 for ti = \xi — . (3.51) 
otherwise 



where 



i=i 



The next step is to decouple a given pair of variables, Xj, Xj, from the other variables in order 
to achieve cancellations between numerator and denominator in the terms on the right side of 
(3.50). As in [13] one defines for i < j < p 

F p ,i = /( min |x p - x fc |), F P:ij = /( min |x p - x fc |). (3.52) 

k<p,kj£i k<p,k<£{i,3} 
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Note that F p j is independent of i and F p ^j is independent of i and j. As in Eqs. (3.13)-(3.15) 
in [13] one has 



77l2 J7i2 77i2 J7i2 ^ JT'2 jp2 T?2 Tp2 ( Q CO\ 



and 



Fj ...F^ > Fj +hj . . . Ff_ hj Ff +l ^ . . . F^ij (3.54) 

N \ I N 

1- (l-/(|x,--Xfc|) 2 ) 1- (l-/(|x,-x fc |) 2 

k=l,k^i,j J \ fc=l,fe^i 

In an analogous way one defines $ Pj j and and obtains the corresponding equations (3.53) 
and (3.54) with / replaced by (p and x by the 2D variable x. 

Consider now the first two terms in (3.50). As in [13] one uses the estimates 

i-l 

f'iU) 2 ^/(Ixi-x,!) 2 , and Fi < /(|xj — Xj|). (3.55) 
j=i 

For fixed i, j one can use (3.53) and (3.54) (and their analogues for $) to separate the contribution 
from the variables Xj and Xj from the rest of the integrand. The part that depends only on the 
other variables cancels between numerator and denominator, while the ij contribution in the 
numerator is 



J (2/'(|x i - Xj |) 2 + v(\ Xi - x,|)/(|x, - x,|) 2 ) 



x Shizifshizjftfdixi - x j \) 2 V GV {x i ) 2 V GV {x j ) 2 dyL i dyi j . (3.56) 



Integration over Z{ and Zj, using Young's inequality, bounds (3.56) by 

W(xi - XjMftxi - Xj\)<p G *{xi) 2 <p G *(xj) 2 . (3.57) 
where W is defined in (3.5). This, in turn, can as in Eqs. (3.18)-(3.20) in [13] be bounded by 

2 J p GP (x) 4 dx,J (3.58) 

with 

J=\j W(x)ip(\x\) 2 dx (3.59) 

given by (3.26) The matching kinetic term, J |V"</?| 2 , given by (3.25) ist derived form the last 
term in (3.50) in the same way. These two terms together give, up to small errors, the correct 
coupling constant, i.e., the factor before J tp GP \x) A dx in the 2D GP functional. The other parts 
of the 2D GP functional as well as the confining energy Nh~ 2 e L follow from the last term in 
(3.50). 

It remains to look at the error terms. These come on the one hand from the ij contribution 
to the denominator, i.e., the last two factors in (3.54) and the corresponding factors with <p 
instead of /. On the other hand they come from the third term in (3.50) and corresponding 



17 



terms with if instead of /. The first mentioned errors have (cf. Eq. (3.21) in [15] and Eq. (3.1) 
in [14]) the form 

^ll^ GP ||Lll^llL(4^ 3 /3), and N \\^ GP W^nb 2 . (3.60) 

Since H^ 013 !!^ ~ L~ 2 and Hs/J^, ~ h~ x , while R = h and b = p~ 1 l 2 \ \n(ph 2 )\~ a , these terms 
are of the order ph 2 and \ln(ph 2 )\~ 2a respectively. For the other errors one again exploits the 
cancellations between numerator and denominator and obtains in the same way as in [15] and 
[14] (cf. Eqs. (3.26) and (3.36) in [13], and (3.4)-(3.5) in [14]) terms of the order {pah) 2 , N' 1 , 
and | m(p/i 2 )| -1 . This can be compared with the estimate for 72-2 i n (3.17) where cancellations 
could not be used. All the mentioned error terms are small in Region II. In combination with 
the bound (3.49) for Region I we can thus state the upper bound in both regions as 

E QM (N,L,h,a) - jV/rV < E%g(N, L, g)(l + o(l)). (3.61) 



4 Scattering length 

As a preparation for the lower bound we consider in this section the perturbative calculation of 
the 2D scattering length of an integrable potential. 

Consider a 2D, rotationally symmetric potential W > of finite range Rq. As discussed in 
Appendix A in [19] the scattering length is determined by minimizing, for R > Rq, the functional 

£rM>]= [ {IWf + iw^l 2 } (4.1) 

J\x\<R 

with boundary condition ip = 1 for \x\ = R. The Euler equation (zero energy scattering equation) 
is 

-AiP + ±Wip = (4.2) 
and for r = \x\ > Ro the minimizer, ip , is 

ipo(r) = ln(r/a sca tt)/ln(i?/a scat t) (4.3) 

with a constant a scat t. This is, by definition, the 2D scattering length for the potential W. An 
equivalent definition follows by computing the energy, 

Er = £r\M = 2vr/ ln(i?/a scatt ) (4.4) 

which means that 

ascatt = Rexp(-2Tr/E R ). (4.5) 

Lemma 4.1 (Scattering length for soft potentials). Assume W{x) = Xw(x) with A > 0, 
w > 0, and w € L 1 (E 2 ), with J w(x)dx = 1. Then, for R > Rq, 

ascatt = R exp I 1 (4.6) 

with 7](X,R) -» for A -» 0. 

Proof. The statement is, by (4.5), equivalent to 

£ R = ±A(l + o(l)) (4.7) 
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where the error term may depend on R. The upper bound is clear by the variational principle, 
taking ip = 1 as a test function. For the lower bound note first that ipo < 1. This follows from the 
variational principle: Since W > the function 4>q{x) = min{l,^o} satisfies £r[V>o] < £rVI>q\- 
Hence the function ipo = 1 — is nonnegative. It satisfies 

-A^ + ±W<p Q = \W (4.8) 

and the Dirichlet boundary condition (p$ = for \x\ = R. 
Integration of (4.2), using that ipo(r) = 1 for r = R, gives 

E R = ±Jwiio = kJw{l-<p ). (4.9) 

Since tpo > we thus need to show that Halloo = o(l). 

By (4.8), and since ipo and W are both nonnegative, we have — Aipo < and hence 

<p (x) < J K (x,x')W(x')dx' (4.10) 

where Kq(x,x') is the (nonnegative) integral kernel of (— A) -1 with Dirichlet boundary condi- 
tions at \x\ = R. The kernel Kq(x,x') is integrable (the singularity is ~ In \x — x'\) and hence, 
if W is bounded, we have ||yo||oo ^ (const. )A||iy|| 00 — O(A). 

If w is not bounded we can, for every e > 0, find a bounded w £ < w with f{w — w £ ) < e. 
Define C £ = \\ 

^ £ ||oo* ^Vithout restriction we can assume that C e is a monotonously decreasing 
function of e and continuous. The function g(e) = e/C £ is then monotonously increasing in e 
(and hence decreasing if e decreases), continuous and g(e) — > if e — > 0. For every (sufficiently 
small) A there is an e(X) = o(l) such that g(e(X)) = A. Then 

llvolloo < (const.)(e(A) + AC e(A) ) = (const. )e(A) = o(l). (4.11) 

□ 

Corollary 4.2 (Scattering length for scaled, soft potentials). Assume Wr : \(x) = 
\R~ 2 w\{x/ R) with w\ > fixed and J w\ = 1. TTten i/te scattering length of Wn t \ is 

a scatt = i?exp(-^±^) (4.12) 

mi/i 77(A) — > for A — > 0, independent of R. 

Proof. This follows from Lemma 4.1 noting that, by scaling, the scattering length of is R 

times the scattering length of \w\. The latter is independent of R. □ 

Remark. If is obtained by averaging a 3D integrable potential v over an interval of length h 
in the z variable, the formula (4.12), together with Eq. (A. 8) in [19], motivates the exponential 
dependence of the effective 2D scattering length (1.17) of v on h/a: The integral A = J W(x)dx 
is h^ 1 J v(x)(i 3 x, which for weak potentials is h~ l Sira to lowest order, by Eq. (A. 8) in [19]. 
Inserting this into (4.12) gives (1.17) (apart from the dependence on the shape function s). This 
heuristics is, of course, only valid for soft potentials v. An essential step in the lower bound in 
the next section is the replacement of v by a soft potential to which this reasoning can indeed 
be applied. 
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5 Lower bound 
5.1 Finite boxes 

Like for the upper bound we consider first the homogeneous case and finite boxes, this time 
with Neumann boundary conditions. The optimal distribution of particles among the boxes is 
determined by using subadditivity and convexity arguments as in [18] and [15]. 

In the treatment of the lower bound there is a natural division line between the case where 
the mean particle distance p -1 / 2 is comparable to or larger than h and the case where it is much 
smaller than h. The first case includes Region II and a part (but not all) of Region I. When 
much smaller than h the boxes have finite extension in the z direction as well. The 
method is then a fairly simple modification of the 3D estimates in [13] (see also Section 4.4 in 
[15]) and will not be discussed further here. 

The derivation of a lower bound for the case that ph 2 < C < oo proceeds by the following 
steps: 

• Use Dyson's Lemma [5, 18] to replace v a by an integrable 3D potential U, retaining part 
of the kinetic energy. 

• Average the potential U at fixed i£l 2 over the z-variable to obtain a 2D potential W. 
Estimate the error in this averaging procedure by using Temple's inequality [37, 29] at 
each fixed x. 

• The result is a 2D many body problem with an integrable interaction potential W which, 
by Corollary 4.2, has the right 2D scattering length to lowest order in a/h, but reduced 
kinetic energy inside the range of the potential. This problem is treated in the same way 
as in [19, 14], introducing a 2D Dyson potential and using perturbation theory, again 
estimating the errors by Temple's inequality. 

• Choose the parameters (size t of box, fraction e of the kinetic energy, range R of potential 
U, as well as the corresponding parameters for the 2D Dyson potential) optimally to 
minimize the errors. 

The first two steps are analogous to the corresponding steps in the proof of the lower bound 
in Theorem 3.1 in [15], cf. Eqs. (3.30)— (3.36) in [15]. It is, however, convenient to define the 
Dyson potential U in a slightly different manner than in [18]. Namely, for R > 2Rq, with Rq 
the range of v a , we define 



The reason is that this potential has a simple scaling with R which is convenient when applying 
Corollary 4.2. Proceeding as in Eqs. (3.30)-(3.38) in [15] we write a general wave function as 




for \R < r < R 



otherwise. 



(5.1) 



n 



*(x!,...,X n ) = /(Xi, 



■n 



) n 



(5.2) 



k=i 



and define F(x±, . . . , x n ) > by 




n 



■n 



)i 2 n s hi z k) 2 dz k . 



(5.3) 



k=i 
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Note that F is normalized if ^ is normalized. The analogue of Eq. (3.38) in [15] is 

n . n 

E / ^^^(i-^IVl^l^min,!^!^)] ]Jdx k 
i=l ^ fe=l 
n n 

+ E / [^l^/l' + ^dxi-x^Dx^^/^i/i 2 ] [] Sfe (z fc ) 2 dx fe , 

-=1 J k=l 

(5.4) 



where Vj denotes the gradient with respect to Xi and <9j = d/dzj. Moreover, xb s is the charac- 
teristic function of the subset B$ C R where s(z) 2 > 5 for 5 > 0, 

a' = a(l - e)(l - 2i?||ds 2 || 0O /(/i<5)), (5.5) 

and /c(i) denotes the index of the nearest neighbor to Xj. When deriving (5.4) the Cauchy 
Schwarz inequality has been used to bound the longitudinal kinetic energy of / in terms of that 
of F, i.e., 

n „ n n „ n 

E / l V !^| 2 II s h{z k ?d* k > E / IVI'^I 2 II dx k . (5.6) 

i=l J k=l i=l k=l 

We now consider, for fixed x±, . . . ,x n , the term 



E / [eldiff + a'URdxi-x^DxBsizkw/hM 2 ] f{s h {z k fdz k . (5.7) 
i=i J k=i 

This can be estimated from below by the expectation value of Ur over Zi at fixed Xj, using 
Temple's inequality to estimate the errors. The result is, by a calculation analogous to Eqs. 
(3.41)-(3.46) in [15], 

_L f n 

(mm-^ > /E[ £ i v f F i 2 + ( 1 - £ )i v ! lF i 2 ^in^i, i -, fc i>«(^) 

i=l 



+\ e w (*i - x i)\ f \ 2 \ n ^ , (5.8) 

j^i k=l 

where W is obtained by averaging ci'Ur over z: 

W(x - x') = 2a" [ Sh {z) 2 s h {z') 2 U R {\x - x'|) X B 5 (z'/h)dzdz' . (5.9) 

JRxR 

Here, a" = a'(l — rf) with an error term r] containing the error estimates from the Temple 
inequality and from ignoring other points than the nearest neighbor to Xj. Moreover, since 
fU R (x)dx = 4tt, E/flflx - x'|) = for |x - x'| > R, and \s(z) 2 - s{z') 2 \ < i?||<9 2 s 2 ||oo for 
\z — z'\ < R we have the estimate 



/ W{x)dx > 
Jm. 2 



^(/ s /(^-f»^l 

8Tia" ||4 R.. 2i 



> — ( \\s\\t-5-j\\d z s 2 \\ O0 ) . (5.10) 
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As we will explain in a moment, the errors, and the replacement of n — 1 by n, require the 
following terms to be small: 

nh?a nR 1 _ R _,. 

!W> IT' £ ' n' * M" (5 - U) 

The rationale behind the first term is as follows. The Temple errors in the averaging procedure 
at fixed xi, . . . , x n produces a factor similar to (2.9), namely 

1 " °' (U) (const.)e//i 2 - (const.)a'(tf) J (5 ' 12) 



with 



„ / n \ m n 

= y [^C/(|x,-x fe(i) |) XB ,(x^ ( . ) /r)j n^(^) 2 ^-, (5.13) 



a=i / j=i 

cf. Eqs. (3.40)-(3.41) in [15]. The analogue of Eq. (3.42) in [15] is 



(U) < (const.)n(n - l)^f , (5.14) 

and (U 2 ) < (const. )nR~^{U) by Schwarz's inequality. Since the denominator in (5.12) must be 
positive, we see in particular that the particle number must satisfy 

eR? 

n(n — 1) < (const.) (5.15) 

ah 

and the error is of the order nh 2 a/eR? as claimed in (5.11). 

The estimate from below on (U), obtained in the same way as Eq. (3.46) in [15], is 

W ^ Y, Ju{\*i-Xj\)XB s {z3/h) (l-E ^-Ixfc-Xil)) f[ Sh (z l ) 2 dz l 

i¥=j \ k,k^i,j J 1=1 

In particular, the second term in (5.11), nR/h, should be small. The requirement that e and 
n _1 are small needs no further comments. 
The potential W can be written as 

W{x) = \R- 2 wi{x/R) , (5.17) 

where w\ is independent of R, with 

y wi(x)dx = 1 (5.18) 

and 

A = J—(l-r/). (5.19) 

Here, r/ is an error term involving 5 and R/(h5) (cf. (5.5) and (5.10)) besides the other terms 
in (5.11). In particular, 5 and R/(h5) should be small. The 2D scattering length of (5.17) can 
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be computed by Corollary 4.2 and has the right form (1.17) to leading order in A. (Recall from 
the remark preceding Eq. (1.19) that R in (4.12) can be replaced by h as long as ca < R < Ch.) 

The Hamiltonian on the right side of Eq. (5.8) can now be treated with the methods of [19]. 
The only difference from the Hamiltonian discussed in that paper is the reduced kinetic energy 
inside the range of the potential W . This implies that A in the error term 77(A) in Corr. 4.2 
should be replaced by X/e, which requires 

4 « 1- (5-20) 
en 

Otherwise the method is the same as in [19]: a slight modification of Dyson's Lemma (lemma 

3.1 in [19]) allows to substitute for W a potential U of larger range, R, to which perturbation 
theory and Temple's inequality can be applied. The modified Dyson lemma is discussed in the 
Appendix. The fraction of the kinetic energy borrowed for the application of Temple's inequality 
as in Eq. (3.16) in [19] will be denoted by e. The errors that have now to be controlled are 

£ ' ^ /£ ' R' eR 2 ln(R 2 /a 2 D )' 

To explain these terms we refer to Eqs. (3.18)-(3.19) in [19] which contain the relevant 
estimates. Substituting e for e and R for R in these inequalities we see first that e and nR 2 /I 2 
should be small. The smallness of R/R guarantees that the "hole" of radius R in the 2D Dyson 
potential has negligible effect. Since the denominator in the Temple error in Eq. (3.19) in [19] 
must be positive we see also that the particle number n in the box should obey the bound 

n(n - 1) < (const.) eln(^ 2 /a| D ), (5.22) 

and the Temple error is bounded by (const. )n£ 2 /{e~R 2 \n(R 2 / a 2 ,^)) . 
We summarize the discussion so far in the following Lemma: 

Lemma 5.1. For all n satisfying (5.15) and (5.22) 

E hox (n,£,h,a) > 2 ™^~ l) \\n{a 2 m /R 2 )\- 1 (l - S (n, £, h, a; e, R, e, R, 6)) (5.23) 

where £ tends to zero together with terms listed in (5.11) and (5.21). 

We note also from Eq. (3.19) in [19] that K(n) = (1 — £ ) is decreasing in n (for the other 
parameters fixed). If K{n) is defined as zero for n not satisfying (5.15) and (5.22) the estimate 
(5.23) hold for all n. 

5.2 Global bound, uniform case 

Using superadditivity of the energy (which follows from W > 0), monotonicity of K(n) and 
convexity of n(n — 1) in the same way as in Eqs. (7)-(12) in [18], one sees that if p = N/L 2 is 
the density in the thermodynamically large box of side length L the optimal choice of n in the 
box of fixed side length £ is n ~ p£ 2 . We thus have to show that it is possible to choose the 
parameters e, R, 5, e and R and £,'m such a way that all the error terms (5.11) and (5.21), as 
well as 5 and R/(hS) are small, while | ln(a2 D /-R 2 )| = | m ( a 2DP)l(l + °(X))- We note that the 
conditions a/h <C 1 and p|ln(/?a2 D )| _1 1/h 2 imply pa 2 ^ — > and hence | ln(pa2 D )| _1 — * 0. 
We make the ansatz 
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and choose £ such that L is a multiple of I with 



P- 1,2 «^P- 1,2 Q^ ■ (5-25) 

Then n = p£ 2 » 1 and R/(h8) = (a/hY~ a ' . The error terms (5.11) are also powers of a/h and 
we have to ensure that all exponents are positive, in particular 

- a' > 0, p - 2 7 > 0, 1 - a - 3/3 - 2 7 > 0. (5.26) 

This is fulfilled, e.g., for 

a = a ' = h p = l 7 = iV (5 - 27) 

with all the exponents (5.26) equal to 1/9. Note also that with this choice Eq. (5.20) is fulfilled. 
Next we write 

e = \ln(pa 2 2D )\- K , R = p' 1 ^ \n(pa 2 2D )\^ . (5.28) 
Then |ln(a^ D /« 2 )| = | ln(a| D p)|(l + o(l)). The error terms (5.21) are 

e=\Hpal D )r, ^ = Q\ph^\Hpat D )\<, ^ = | ln(pa 2 D )|^, (5.29) 



and 



nt 2 ( a \ ~ 4 7 



7|hl( "°- ,r(1 """ 2£) (l + 0(.n|ln(p„I D )|)). ,5.30) 
Since (a//i) -47 | ln(/9a 2 D )| _47 = 0(1), the error term (5.30) can also be written as 

nf 2 

, 521 L, 2 : = 0(1)1 ln^r^-^Hl + 0(ln | ln(^ D )|)). (5.31) 
eR 2 ln(R 2 /<%,) 

The condition p/i 2 < C is used to bound R/R in (5.29). Namely, 

(p/i 2 ) 1 / 2 )!^/^)^ < (const.) {h/af, (5.32) 



so 



We choose now 



Then 



C=^, « = (5.34) 

/3-C = ^ and 1- K -2C-4 7 = 1 (5.35) 

This completes our discussion of the lower bound of the energy in the thermodynamic limit, 
i.e., in Theorem 1.4, for the case pti 2 < C < oo. As already mentioned, the case ph 2 3> 1 can 
be treated with the 3D methods of [18, 13], taking care to retain uniformity in the parameter 
h/l. The analogous problem for the reduction form 3D to ID is discussed in detail in [15], cf. in 
particular the lower bound in Theorem 3.2 in [15]. Since no new aspects arise in the reduction 
from 3D to 2D we refrain from discussing the case ph 2 3> 1 further here. 
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5.3 Global bound in a trap and BEC 

We now discuss the case when the system is inhomogeneous in the x variables, i.e., the lower 
bound for Theorem 1.1. A possible approach, analogous to that of [13], is to modify Eq. (5.2) 
and write a general wave function as 



*(xi, . . . ,x n ) = /(xi, ,x n ) Yl f GP (x k )s h (z k ) 



(5.36) 



k=i 



This is always possible because ip and Sh are strictly positive. As in [13] the task then becomes 
to minimize a quadratic form in / involving the GP density p GP (xj) = ip (xi) 2 both as a weight 
factor in the integration over the Xj's and also as a replacement for the external potential Vl(xj). 
There is, however, a neat alternative approach to deal with inhomogeneities due to R. Seiringer 
[36] that is somewhat simpler and that we shall follow with appropriate modifications. This 
approach leads also quickly to a proof of Theorem 1.3. 

5.3.1 GP case 

We start by writing again the wave function in the form (5.2) but this time with n = N. With 
H including the trapping potential Vl, cf. (1.3), and using the symmetry of F and / we obtain 
as in Eq. (5.4) 

N 1 f N 

(*\H\V) >T + N V L ( Xl )F 2 J] dx k 

J k=i 

N 

+N / [eld^ + a'UR^-^DxB^m/hM 2 ] l[s h (z k ) 2 dx k , 
J i — i 



where 



T = N 



r N N 

/ V ; VI/-1J ( / X/ • / Vf\|/| 2 J] dx fc 

■/|xi-x fc( i)|<fl k=1 J\xi-Xk(i)\>R k=l 



(5.37) 



(5.38) 



\ x i- x k(i)\<R 

Note that we have here not yet used Eq. (5.6) to bound the parallel kinetic energy of / in terms 
of that of F. The reason is that we want to prove BEC for the wave function ^ and not only 
for F. 

The next step is to split the kinetic energy T, for given R> R and e > 0, into two parts with 
the main contributions coming from length scales > R and < R respectively in the longitudinal 
directions. To write this in a compact way we introduce the notation 



X = (x2, ■ ■ ■ , xn), dX = dx2 ■ ■ ■ <1xn', Z = (z±, zn), dZ = dz\ ■ • • dz^ 



(5.39) 



and 



R Ft 

x — {xi : \xi - x ki i) | < R}, tl x ' = {xi : R < \xi - x fe(1) [ < R}, 



The splitting of the kinetic energy is 



T > T > + T < 



.v - {xi : \xi -x fc( i)| > R}. (5.40) 

(5.41) 
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with 



dXdZ, (5.42) 



dX. (5.43) 



t> = n[ e l\ \vh\ 2 + ^ [ iv 1 , 1 ^! 2 + (i - £) f \vh\ 

J R tN-2 [ 2 J n <R 1 1 1 2 y n «.A 1 1 1 ^ 2^>*' 1 1 

t< = n[ £ h-i) f |v)f| 2 + (i - J) / iv)fi 2 + / J|vJf| 

Here we have used (5.6) for T < , but not for T > . We then add and subtract 
8nN 2 g J R2N p GP (xi)|F| 2 from the right side of (5.37) and, using (5.41), estimate it by E > + E < 
with 

E> = T> +N j V L {xi)F 2 + 8vriV 2 5 J p GP (xi)|F| 2 , (5.44) 



N 



E <=T< + N f [el^/p + a'l/fldxi-Xfc^DxB,^!)//!)!/! 2 ] f[ *h(*k) 2 <kfc 

fc=i 

-8vriV 2 5 y p GP (xi)|F| 2 . (5.45) 

The two terms are now considered separately. The first remark is that ip GP (x) is the ground 
state wave function of the 1-body Hamiltonian 

-A 11 + V L (x) + 8nNgp GP (x) (5.46) 

and the energy is equal to the GP chemical potential 

p GP (Ng) = E GP (l,L,Ng)+4iTNg J <p GP (x) A dx, (5.47) 

cf. Section 2 in [14]. Hence, if T > were replaced by the full kinetic energy, J |v|\l/| 2 , then E > 
would be bounded below by Np, GP (Ng) = E G ^(N,L,g) + ATrN 2 g J ip GP (x) 4 dx. Using a variant 
of the generalized Poincare inequalities of [16] it is shown in [35, 36] that this estimate also 
holds for E > up to small errors. Even a stronger result is true (cf. [36], Eq. (51)), because a 
contribution form the one particle density matrix in the lontitudinal direction, 



l\{x,x')= I jq,((x,z),(x',z))dz 

can be included: Provided R 2 N —>■ for N — > oo, and Ng is fixed, 

E>/N > (» GP (Ng) + CTr [^ 7 |(1 - \ip GP )(ip GP \)] ) (1 - o(l)) 



(5.48) 



(5.49) 



where the term o(l) goes to zero if e, e — > go to zero after the limit N — > oo has been taken. 

The term E < is treated in the same manner as in the homogeneous case by introducing 
boxes, A Q , of side length t in the longitudinal directions. Defining 



p a = sup p GP (x) 



we have (cf. Eq. (60) in [36]) 



E< > inf y^(E(n a ,e) - 8itNgp a n a ) 



a 



(5.50) 



(5.51) 
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with E(n,£) the right side of (5.23). (It depends also on the parameters h,a;e,R,e,R,8.) The 
infimum is taken over all distributions of the N particles among the boxes and n a denotes the 
particle number in box a. Choosing the box length as in Eq. (5.25) with p replaced by the mean 
density p ~ N/L, and arguing exactly as in [36], Eqs. (71)-(78) we obtain 

E< > -inNg J \ip GP (x)\ 4 dx(l + o(l)). (5.52) 

Adding (5.52) and (5.49) we obtain the bound 

- ^ > [e gf (1,L,N 9 ) + CTr[i 7 |(l - |^ GP )(^ GP |)] (1 - o(l)). (5.53) 

If ^ is the ground state wave function ^o, then (5.53) together with the upper bound (3.61) 
proves both the convergence of the energy as claimed in Theorem 1.1 and also the convergence of 
the density matrix j^(x, x') to cp GP (x)<p GP (x') (in trace norm) as N — > oo for the case that Ng is 
fixed (or at least bounded). The convergence of the full density matrix 7$ (x, x'), scaled by h in 
the variables z and z', to ip GP (x)tp (x')s(z)s(z') follows from a simple energetic consideration: 
The energy gap between the lowest and the first excited level of —d% + is large compared to 
(^f\H\^f)/N — j-? in the limit considered. Hence the contributions to the density matrix form 
the excited states in the perpendicular direction must vanish. 

Remark. The above proof of BEC can be extended to allow a slow increase of Ng to infinity 
with N, but this improvement is restricted by the available estimates of the error terms in (5.49) 
and (5.53) and is only marginal. 

5.3.2 TF case 

To complete the proof of Theorem 1.1 one must also consider the case that Ng — > oo, i.e. the 
'Thomas- Fermi' limit. Here the starting point is again (5.37), this time with T replaced by the 
smaller quantity 



T 



f N N 

£ / |V|F| 2 l[dx k + / \vlF\ 2 ]Jdx k 

J\xt—Xi.ri<i\<R ;„_1 J Ixi — It/il \>R , -i 



(5.54) 



I |xi-x fc (i)|<R k=1 J \xi-x k{1) \>R k=1 

cf. (5.6). We then use the TF equation 

pTF{x) = 8^Ng- [flTF - VL{x)]+ (5 ' 55) 

to trade the potential Vl(x) for the TF density p TF (x) which is the square of the minimizer 
of (1.10) with the kinetic term omitted. The TF chemical potential p is determined by the 
normalization condition J p TF = 1 and satisfies in analogy to (5.47) 

p TF (Ng) = Ej p (l, L, Ng) + iirNg f V TF (x) 4 dx, (5.56) 

From (5.55) it follows that 

V L {x) > p TF - 8Tr(Ng)p TF (x). (5.57) 

Both p TF and p TF depend on Ng and scale in a simple way if the potential Vl{x) is a homoge- 
neous function of some degree p > 0, cf. Eqs. (2.20)-(2.21) in [15]. We can now apply the box 
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method to the minimization of 

N 



T + M TF - 8n(Ng)N f p TF ( Xl )F 2 J] dx k 
^ k=i 

N 

+ N / [el^/l 2 + a'C/fldx! - ^ m \)xB s {z k{l) /h)\f\ 2 ] \[ s h (z k ) 2 dx k , (5.58) 
J k=i 

obtaining in analogy to (5.51) and using (5.56), 

(5.58) > E%g(N,L,g)+4ir(Ng)N [ p TF (x) 2 dx + inf V(£(n Q , £) - &irNgpl F n a ). (5.59) 

J n a ^ 

The remaining steps towards the estimate 

(5.58) > Ej F (N,L,g)(l - o(l)) (5.60) 
are now exactly like in Section 4 in [14]. 

6 Conclusions 

We have investigated the dimensional reduction of a trapped, interacting Bose gas when the trap 
potential is strongly confining in one direction. Starting from the many-body Hamiltonian with 
repulsive, short range interactions we have shown rigorously how an effective 2D description 
of the ground state energy and density emerges and how the parameters of the 2D gas relate 
to those of the 3D gas. Two parameter regimes can be distinguished: One where the gas 
retains some of its 3D character despite the tight confinement and another where the situation 
is manifestly two-dimensional with a logarithmic dependence of the coupling parameter on the 
density. Moreover, we have shown that the trapped gas is Bose-Einstein condensed in the ground 
state provided the coupling parameter Ng stays bounded. 

A Appendix: A modification of Dyson's lemma 

Let W(x) > be a rotationally symmetric potential in 2D with W(x) = for |x| > R. Denote 
by Br the ball of radius R around the origin. For < e and R' > R define 



E R ' tE = min 



f [e\V<P{x)\ 2 + \W{x)\<p(x)\ 2 } dx+ j \V4>{x)\ 2 dx (A.l) 

J B R J B R/ \B R 



where the minimum is taken over 4> G H l {B R ,) with 4>(x) = 1 for |x| = R' . 
Lemma A.l. ^ 

ER ' £ = HR'/R)l 2n/E R>£ (A ' 2) 

Proof. For any c > the minimum of the first integral in (A.l) with boundary condition 4>{x) = c 
for \x\ = R is c 2 E R e and the minimum of the second integral with the same boundary condition 
at \x\ = R and <p(x) = 1 for |x| = R 1 is 2ir \n(R' / R) / (\w.(R/b)) 2 where b is determined by 
ln(i?/6)/ ln(R'/b) = c. Adding the two contributions and minimizing over c gives (A. 2). □ 
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Lemma A. 2 (Modified Dyson Lemma). Let W be as above, let R > R and let U(r) > be 
any function with support in [R, R] satisfying 

r R 

2tt / U{R')E R } £ R'dR' < 1 (A.3) 

J R 

with Eft £ as in (A. 2). Let B C R 2 be star-shaped with respect to 0. Then, for all functions 
[ e\V<p(x)\ 2 + \w{x)\<f){x)\ 2 dx+ [ \V<p{x)\ 2 dx 

[ f/(|x|)|0(x)| 2 dx. (A.4) 
Jb 



lBn(B A \B R ) 

> 



Proof. The proof is very similar to that of Lemma 3.1 in [19]. In polar coordinates, r, 6, one has 
|V0| 2 > \d(p/dr\ 2 so it suffices to prove the analogue of (A.4) for each angle 6 € [0, 2ir). Denote 
(j)(r, 6) simply by f(r), and let R(9) denote the distance of the origin to the boundary of B along 
the ray 0. It suffices to consider the case R < R{9) (here, W > is used) and the estimate to 
prove is 



rR rmin{R,R{6)} 

I {e\df(r)/dr\ 2 + \W {r)\f {r)\ 2 } rdr + / \df(r)/dr\ 2 r dr 

Jo ~ JR 

min{R,R(8)} 



> / U(r)\f(r)\ 2 rdr. (A.5) 

JR 

For the given value of 9, consider the disc B R ^ centered at the origin in M 2 and of radius R(9). 
Our function / defines a rotationally symmetric function, x i— > /(|x|) on B R ^, and (A.5) is 
equivalent to 



f \e\Vf(\x\)\ 2 + ^(r)|/(|x|)| 2 | dx + [ |V/(|x|)| 2 dx 

> j U(\x\)\f(\x\)\ 2 dx (A.6) 



B R(0) 



If R < R' < min{i?, R(9)} the left side of (A.6) is not smaller than the same quantity with 
Br(6) replaced by the smaller disc Br>. (Again, W > is used.) According to (A.l) this integral 
over Bri is at least Er> >£ \f(R')\ 2 . Hence, for every such R', 

rR rmm{R,R(6)} 

/ {e\df{r)/dr\ 2 + \W{r)\f{r)\ 2 } rdr + / \df(r)/dr\ 2 r dr 

Jo Jr 

> (2vr)- 1 ^, £ |/( J R')| 2 (A.7) 

The proof is completed by multiplying both sides of (A.7) by 2ttU(R')E^} £ R' and, finally, 
integrating with respect to R' from R to min{^, R{9)}. □ 

A convenient choice for the Dyson potential U is 

m = [ XR<r<R (Ag) 

10 otherwise , 
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with 

v{R) = 2ir / E RI>£ R' dR'. (A.9) 
Jr 

In the case considered in this paper the potential W is integrable with coupling constant 
A ~ a/h(l + o(l)) while e ~ (a/Zi) 1 / 9 . In particular X/e — > as a/h — > 0. From Section 4 and 
the discussion of Eqs. (5.19)-(5.20) we thus conclude that in our case 

E R>£ = 2vr/ ln(i?/a 2D )(l + o(l)). (A.10) 

By Eq. (A. 2) we then also have 

£^, e = 2vr/ln(i?7a 2D )(l + o(l)). (A.ll) 

and thus 

= p 2 ln(i? 2 /«2D) + (A.12) 

Thus, for R and e chosen as in (5.24) and (5.27), the modified Dyson Lemma gives, up to 
negligible errors, exactly the same result as the standard 2D Dyson Lemma for a potential with 
scattering length a 2 D- 
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